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Preface 



This report is a collection of material that has been used 
in courses on search, detection and localization modeling. Its 
organization follows to some extent material by S. M. Pollock in 
Selected Methods and Models in Military Operations Research which 
is listed in the report bibliography. The report is not intended 
to be a text on these subjects. In particular, in some areas it 
does not provide the depth of coverage that is found in the book 
Search and Detection by Alan R. Washburn which is cited in this 
report as Reference 22. 

In the second revision, some typographical and other errors 
have been corrected and some changes have been made to several 
sections of the report. 
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I. Detection Models and Signal Detection Theory 

Signal detection theory is the basis for analyzing the 
detection models that are described in this report. In signal 
detection theory, the decision making portion of a detection 
system is called the receiver and a detection experiment is the 
observation by a receiver of input data generated during some 
time interval. The data that is related to a target is called 
signal. The data that is not related to the target is called 
noise. In general, the target data is associated with a 
localization region that in some cases is called a resolution 
cell. When a detection experiment is performed, either the event 
Hq = (the receiver input is noise} or its complement 

= {the receiver input is signal and noise} will occur. In 
the first detection models that are described here, after 
analysis of the input data by a receiver, either the event 
Dq ={the receiver decides the input is noise} or its complement 
D^ = (the receiver decides the input is signal and noise} will 
occur. Detection models for which D^ is the complement of Dq 
are called binary detection models or forced choice detection 
models. Four events which are important in binary detection 
models are indicated in the Venn diagram of Figure 1. 

The Venn diagram emphasizes a decision problem that is 
associated with a receiver that can be modeled using a binary 
detection model. The problem is this: Under what conditions 

should the event occur? That is, under what conditions 

should a receiver decide that the input data indicates a target 
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was present in an observed localization region either at the time 
or prior to the time of the observation? 



Ho Hi 



Dq ^ Hq 


Dq n Hi 


Di n Hq 


Di n Hi 



Figure 1. Four events of importance in binary detection models. 

In the detection model descriptions that follow, the 
following notation and terminology are used: pf = P(Di|Ho), the 

probability of Di given Hq, is called the false alarm 
probability; p^j = P(Di|Hi), the probability of Di given Hi, 
is called the detection probability and P = P(Hi), the 
probability of Hi, is called the prior probability. It is the 
probability that a target will be in the localization region when 
the input data is generated. 

In the detection models, the input to a receiver is 
determined by a stochastic process that has the following 
characteristics: It is a random noise process when there is no 

target data and it is a random noise process plus a signal 
process when there is target data. Although the receiver input 
process in some cases might appear to be determined by a 
continuous parameter stochastic process, because of the finite 



2 



amount of information (unique data) contained in a bounded 
sequence of finite length, a discrete parameter stochastic 
process is sufficient to determine the receiver input in these 
cases. This is established formally by the stochastic sampling 
theorem. Consequently, in these models, the input to a receiver 
is determined by a sequence of random variables Y^, ♦♦♦ ,Yj„ and 
an observation yields a sequence of values yi/***Ym' 

Three detection models are described in Section III. In the 
first model, the signal process is a deterministic process. In 
the second and third, the signal process is a random process. To 
define a random noise process or a random signal process, only 
the joint distribution of the finite sequence of random variables 
that determine the process needs to be specified. If the signal 
process is a deterministic process, the signal values can be 
determined before an observation is performed. To define the 
process in this case, only these values need to be specified. 



II. Decision Criteria 

To simplify the discussion of decision criteria and 
decision rules, a receiver's input will be assumed to be 
determined by a single decision random variable Y. In this 
case, the input process in determined by the conditional 
distribution function Fy(y|Ho) when the input is noise alone 
and by the conditional distribution function FyCylH^) when the 
input is signal plus noise. 

The condition that a receiver's input is recpaired to 
satisfy in order that the event will occur can be specified 

in terms of a decision rule. For the assumed case, a decision 
rule is a rule which determines for every observable value of Y 
the decision that the receiver is to make. The decision rule can 
be considered to be a function 4>(y) which relates each 
observable value y to one or the other of the following two 
commands: do = "decide that the receiver input was noise" and 

d]^ = "decide that the receiver input was signal and noise". 
Choosing a decision rule <t>{y) defines a set n of observable 
values of Y such that the event = { Y e n } . 

The problem which was considered in Section I can now be 
restated in the following way; What criterion should be adopted 
in order to determine a decision rule or, equivalently, its 
corresponding set D ? A desirable characteristic for a criterion 
is suggested by the following argument: Consider the odds in 

favor of given y is observed. That is, consider 
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P(H]^|y = y)/P(Hq|y = y) . One might expect that y would be a 
member of the set n if and only if y made this ratio equal to 
or greater than some value k. But this is equivalent to defining 
n as follows; n = { y : L(y) > K } where L(y) is the 
likelihood ratio associated with an observed value y and K is 
a constant related to the constant k. This suggests that 
choosing an optimum criterion is equivalent to choosing an 
optimum value for K. Four specific decision criteria are defined 
next in terms of K. For each criterion, fl has the above 
form. But for each criterion the choice of K is different. 

The decision criteria are: 

1. The Neyman-Pearson Criterion: Choose n so that p^ 

is a maximum subject to the constraint that pf < a where a is 
a specified value. For a continuous decision random variable, 
the constant K is chosen so that pj = a. 

2. The Bayes Criterion: Choose fl so that the expected 

cost of a receiver's decision is a minimum. For a continuous 
decision random variable, K = [ (c^o'^^oo)/ ('^Ol'^ll) 3 ’ 

^10 ^ ^00 *^01 ^ *^11 where cj^j is the cost of Dj^ n Hj . 

3. The Ideal Observer Criterion: Choose n so that the 

probability that the receiver makes an incorrect decision is a 
minimum. K = (1-P)/P for a continuous decision random variable. 

4. The Minimax Criterion: Choose O when P is unknown so 

that the maximum expected cost of a receiver's decision is a 
minimum. If C]^q > ^00 <^01 ^ ‘^ll' then 
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K = [ (C]^o“*^Oo) / (^01“^ll) 3 * ( ^”^*) /^* ^ continuous decision 

random variable. Here, P* is the value of the prior probability 
P that would make the expected cost of a receiver's decision a 
maximum if P were known and the Bayes Criterion were used. 

If a model which specifies the conditional distributions 
FY(yl^o) FY(yl^l) ^ decision rule are adopted, then 

the value of pf and the value of p^j are determined. This 
pair of values (PffPd) called a receiver operating point. 

If the decision rule results from using a likelihood ratio 
criterion such as one of the four listed above, then it will 
involve the parameter K since n = {y: L(y) > K}. And, for a 
given value of K, since n uniquely determines the pair 
(Pf /Pd) / ^ single operating point results. By varying K, a set 
of operating points can be generated which determines a receiver 
operating characteristic curve or ROC curve. Different ROC 
curves can be produced by changing either one or both of the 
conditional distributions which implies either a change in the 
signal process or a change in the noise process. 

A decision rule which results from using a likelihood ratio 
criterion in a model in which the input process is determined by 
a set of m random variables can be expressed in terms of a set 
n as follows: n = { (y^, ••• ,yj^) : L(yi, ••• ,yjj,) > K } where 
K is specified in the same way that it is when m = 1. 
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Ill . Three Binary Detection Models 

Three detection models are examined in this section. 

For the first two detection models, the input stochastic process 
for an observation is defined by a time sequence of continuous 
random variables. The random variables represent a sample from a 
continuous parameter stochastic process which is sampled at times 
such that the random variables are independent. For the third 
detection model, the input stochastic process is a counting 
process and it is defined by a single discrete random variable 
that is equal to the number of events that are counted during the 
observation. 

Model I: In the first detection model, a sampled noise 

value is a value of a normally distributed random variable with 
mean zero and with known variance a*. And a sampled signal 
value is a known value of a deterministic variable. Thus, the 
input process corresponding to an observation consists of some 
number m of independent normal random variables ••• ,Yjjj 

each with variance a*. And, for i = 1,2, ••• ,m, when a signal 

is not present the mean of Yj^ is zero and when a signal is 
present the mean is s^. The result of using a likelihood ratio 
decision rule in the model can be expressed in terms of a random 
variable Z. This random variable is called a crosscorrelation 
statistic and it is defined by Z = S sj^'Y^ where the sum index 
1 = 1,2, ••• ,m here and in the remainder of this section. 

However, it is more convenient to express the result in terms of 
a statistic V which is defined by V = In terms of this 
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random variable, the conditional probabilities pf and p^j are 
given by; pf = 1 - $(v*) and p^j = 1 - #(v* - d^) where § 
symbolizes the standard normal cumulative distribution function, 

V* = (l/Cg) • (a*ln K + (1/2) 2 s|^) , is determined by the decision 

2 1 , , , 
rule and d = 2 Sj^/a is called the detection index. 

Often, the input stochastic process represents a quantity 

whose square is proportional to power. In such a case, the 

. . . 2 
average receiver input power is the random variable 2 Yj^/m. If 

a signal is not present, the expected average receiver input 

power is N = 2 a*/m = a* where N is called the noise. The 

. . . . * 
average receiver input signal power is S = 2 sj^/m where S is 

called the signal. In these terms, d = m* (S/N) where S/N is 

called the signal-to-noise ratio. 

If a receiver's input can be considered to be a time 

sequence of continuous voltage values such as in the case of a 

sonar receiver, in some cases a frequency representation can be 

used that involves the concept of receiver bandwidth. In these 

cases, the noise process is assumed to be such that m = t/5t 

where t is the integration time (the duration of an 

observation) and 5t is the time between samples with 

St = 1/[2(BW)] where BW is the bandwidth and St is 

determined by the sampling theorem. This implies that the 

detection index can be written as d = 2t* (BW) • (S/N) . By 

defining Ng as the power spectral density where Ng = N/BW^ the 

detection index can also be written as d = 2t* (S/Ng) . 
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In Reference 2, the conditions required by the first model 
are called Case I and in the following sections the first model 
is called the Case I model. A receiver that processes data 
such that it would implement a likelihood ratio decision rule 
under the conditions of the first model is called a matched 
filter or crosscorrelation detector. If the description of the 
input noise is adequate, a Case I model can be used to obtain an 
estimate of an upper bound on a detection system's performance, 
since all the information necessary to define the signal is 
assumed to be known. 

Model II; In the second detection model, a sampled noise 
value is a value of an independent normal random variable with 
mean zero and with known variance a*. And a sampled signal 
value is an independent random variable with mean zero and with 
known variance a|. Thus, the input process corresponding to an 
observation consists of some number m of independent normal 
random 

variables each with mean zero and each with variance 

a* when a signal is not present and with variance a* + o| when 
a signal is present. The result of applying a likelihood ratio 
decision rule in this model can be expressed in terms of a 
statistic X which is defined by X = S Y^. When a signal is 
not present, the statistic X/N has a chi-square distribution 
with m degrees of freedom. When a signal is present, the 
statistic X/(N+S) has a chi-square distribution with m 
degrees of freedom. So, in terms of these statistics, the 
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conditional probabilities pf and p^j are defined as follows: 
Pf= P(X^ > x*/N) and p^ = P{X^ > (x*/N) • [ 1/ ( 1+S/N) ] } 
where X^ is a chi-square random variable with m degrees of 
freedom, x* is a number which is determined by the decision rule 
and S/N is the signal-to-noise ratio. A receiver that would 
implement a likelihood ratio decision rule under the conditions 
of the second model is called an energy detector or square law 
detector. 

The mean of a chi-square random variable with m degrees of 
freedom is m and the variance is 2m. By the central limit 
theorem, as the number of degrees of freedom of a chi-square 
random variable becomes large it can be approximated by a normal 
random variable with the same mean and variance. Hence, for a 
sufficiently large sample size m or equivalently, for a 
sufficiently large time bandwidth product t* (BW) , the 
conditional probabilities pf and p^j can be approximated by 
Pf = 1 - f(v*) and Pd = 1 - §{ [1/ (1+S/N) ]• (V* - d^)) where the 
threshold value v* = (x* - mo* ) / [ (2m) * ] and d = t* (BW) (S/N) * 
if the concept of bandwidth is applicable. If, in addition, the 
noise N is significantly larger than the signal S, that is, if 
t* (BW) » 1 and S/N « 1, then pf and p^j can be 
approximated by: Pf = 1 - #(v*) and p^j = 1 - *(v* - d^). 

In Reference 2, the conditions required for this 
approximation are called Case II and in the following sections 
the limiting form of the second model is called the Case II 
model . 
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Model III: In the third detection model, a sampled noise 

value and a sampled signal value are values of independent random 
variables that are determined by independent Poisson processes 
that is observed for a time interval t. The noise process is 
characterized by a counting rate a, the signal process is 
characterized by a counting rate Ug and the noise and signal 
processes are additive. This implies that when the input is 
noise alone, the input is a Poisson random variable with 
parameter a*t and when the input is signal and noise, the input 
is a Poisson random variable with parameter (a + ag)*t. 

For a likelihood ratio decision rule, pf = 1 - P(y*;a*t) 
and P (3 = 1 - P[y*;(ct + ag)‘t] where y* is a threshold value 
that is determined by the decision rule and P(y;9) represents 
the Poisson cumulative distribution function with parameter 0. 

When 9 is large, the cumulative distribution function can 
be approximated by the cumulative distribution function of a 
normal random variable that has the same mean and variance. 

Using this approximation for cases where a*t is large, since 
both the mean and variance of a Poisson random variable are equal 
to 9, Pf = 1 - #(v*) and p^j = 1 - #{[1/(1 + ag/a)^](v* - d^)} 
where v* = (y* - at)/(at)5 and d = a*t*(ag/a)*. In addition, 
if the signal counting rate is significantly smaller than the 
noise counting rate, that is, if a*t » 1 and oLq/cl « 1, then 
Pf and P(j can be approximated as follows: pf = 1 - #(v*) and 

Pd = 1 - #(v* - dM . 
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The third detection model might be used to describe a 
receiver whose input for an observation is the number of photons 
counted by a radiation detector in situations where a*t, the 
expected number of counts when no signal is present, is of the 
order of thirty or more. 

When a likelihood ratio decision rule is used in the three 
models discussed above, for the first model and under limiting 
conditions for the second and third models, the following result 
is obtained; pf = 1 - #(v*) and p^j = 1 - #(v* - d^) where the 
definition of v* depends on the noise power N for the first 
and second models. For a sonar receiver described by the first 
model, that is, by the Case I model; d = 2t* (BW) (S/N) . For a 
sonar receiver described under the limiting conditions for the 
second model, that is, by the Case II model, d = t • (BW) (S/N) * . 
So, in either a Case I model or a Case II model of a sonar 
receiver, the detection index d is a function of the time 
bandwidth product t* (BW) and the signal-to-noise ratio S/N. 
Since sonar equations relate S/N to system, target and 
environmental parameters, a sonar equation can be used to relate 
S/N to these parameters in a model of a sonar receiver. 
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IV. General Detection Models 



The detection models that have been considered to this 
point are based on binary detection theory. After each 
observation, a receiver decides either that the input 
corresponding to the observation was noise or else it decides it 
was signal plus noise. However, in some detection systems this 
decision is delayed. In a computational sense, a model of such a 
detection system is generally complex relative to a binary 
detection model. To illustrate this, consider an active sonar 
system whose receiver includes an operator. Suppose the 
probability that the operator will detect a target echo has been 
determined in a laboratory experiment in which the operator was 
required to decide after each input corresponding to a resolution 
cell that either the input was a target echo (signal) and noise 
or the input was noise alone. In addition, suppose that under 
operational conditions the operator normally delays this 
decision. Then, in general, the probability that the operator 
will decide that the input corresponding to a resolution cell 
that contains a target is a target echo and noise will not be 
equal to the probability of the event in the forced choice 
experiment. And, in addition, the probability that the operator 
will decide the input corresponding to a resolution cell that 
does not contain a target is a target echo and noise will not be 
equal to the probability of this event in the forced choice 
experiment. Consequently, in general, the value of both p^j and 
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Pf for an operational environment will be different than that 
for the laboratory environment. 

One model that has been proposed to deal with this kind of 
situation defines the event that a receiver decides that the 
input corresponding to a resolution cell is signal and noise to 
be equivalent to the event that out of n consecutive 
observations at least k of them would result in the decision 
that the input was signal and noise in a forced choice 
experiment. The model is said to be based on an k-out-of-n 
detection criterion. With this criterion, the probability that a 
target will be first detected on the observation can be 

found as follows: Determine the 2^ sequences of forced choice 

responses that could result for a sequence of j consecutive 
observations. Next, determine the probability of occurrence for 
each sequence that first satisfies the k-out-of-n detection 
criterion on the observation. The probability of first 

detection on the observation is equal to the sum of these 

probabilities. The cumulative probability of detection at the 
jth observation is the sum of the probabilities of first 
detection on the i^^ observation for i=l,2, ••• ,j. 
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V. Signal-to-Noise Ratio Detection Models 

In some radar and sonar detection models, for a 
specified value of pf, a minimum acceptable value of p^ is 
defined. This minimum acceptable value of p^ and the specified 
value of Pf define what can be called a minimum acceptable 
signal-to-noise ratio (S/N)jj, if p^j is a nondecreasing 
function of signal-to-noise ratio. In some sonar detection 
models, (S/N)jj, in decibels is called the detection threshold 
DT. In symbols, DT = 10 log(S/N)jjj. If the minimum acceptable 
value of P(j is .5, then DT is usually called the recognition 
differential RD. The difference between the signal-to-noise 
ratio in decibels and RD (or DT) is called the signal excess 
SE. In symbols, SE = 10 log(S/N) - RD. 

One interpretation of signal excess is that for a 
localization region containing a target detection occurs with 
probability one if SE > 0 and with probability zero if SE < 0. 
This interpretation provides the basis for defining detection in 
the three encounter detection models that are discussed in 
Section VII. A more consistent interpretation of signal excess 
is: If SE > 0, then the probability of detection p^j is 

greater than or equal to the minimum acceptable value (.5 if 
recognition differential RD is used to define signal excess) . 
For cases where p^j increases rapidly with signal excess in the 
neighborhood of zero signal excess, the two interpretations may 
be operationally equivalent. For a discussion of this point as 
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well as a discussion of an operational case in which receiver 
decisions are delayed, see Reference 3. 

Signal excess (signal-to-noise ratio) detection models 
provide a basis for general detection models, in particular, 
models that describe nonstationary noise and signal processes and 
randomly changing decision rules. This is illustrated by the 
models described in Section VII. In addition, signal excess 
models provide a basis for delayed receiver decision models. 

This is illustrated by the active sonar detection models in both 
Reference 4 and Reference 5 that are based on a k-out-of-n 
detection criterion. In all of these models, the signal-to-noise 
ratio and the recognition differential are random variables. 

Using X(t) to represent a random variable corresponding to 
an index time t and a subscript to identify the random variable 
in such models, for a passive sonar receiver, the signal-to-noise 
ratio in decibels associated with a decision at the index time 
is: XsL(t) - X^L(t) - [XNL(t) - XDi(t)]. In this expression, 

SL represents source level, TL represents transmission loss, 

NL represents noise level and DI represents directivity index. 
Since signal excess SE is defined to be the difference in 
decibels between the signal-to-noise ratio and the recognition 
differential (or detection threshold) , it too is a random 
variable and, for any decision time t, one can write: 

(1) XsE(t) = XsL(t) - XTL(t) - [XNL(t) - XpjCt)] - XRo(t). 

The distributions of the random variables on the right side of 
Equation 1 determine the distribution of the signal excess. In 
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the passive sonar detection model described in Reference 6, 
f and, in effect, XjjL(t) are normally 

distributed random variables while X-pL(t) is a uniformly 
distributed random variable. In the three signal excess models 
that are described in Section VII, all of the random variables in 
Equation 1 are normally distributed. 

It is sometimes convenient to write Equation 1 as follows: 

(2) XsE(t) = SE(t) + X(t) . 

In Equation 2, SE(t) is the expected value of the signal excess 
determined by the following expected value equation: 

(3) SE(t) = SL(t) - TL(t) - [NL(t) - DI(t)] - RD(t) 

where each term on the right represents the expected value of the 
indicated random variable and X(t) is a random variable that 
determines the stochastic character of the signal excess. Since 
SE(t) is the mean of XsgCt); by Equation 2, the mean of X(t) 
is equal to zero and the standard deviation of X(t) is equal to 
the standard deviation of XgE(t). If a represents the 
standard deviation of Xgg(t) and the random variables on the 
right side of Equation 1 are statistically independent, then 

~ “^SL ^NL ^DI *^RD' This relation has been used to 

determine a standard deviation for the signal excess in 
operational models. 
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VI. General Encounter Models 



A basic problem associated with search modeling is that 
of determining the probability that a target will be detected by 
a detection system during an encounter with one or more detection 
systems. In the encounter models that are considered in this 
report, during a search, observations are made of a series of 
localization regions. The probability of detection on an 
observation is n . The probability of a false alarm on 

an observation is P(Di n Hq) • In these models, the time to 
resolve a false alarm is ignored. However, p^j and Pf are 
assumed to be determined by some criterion such that Pf is an 
operationally reasonable value. 

Using the order number of a decision rather than its time as 
an index relative to detection decisions for localization regions 
that contains a target and a random variable N to represent the 
decision order number at which detection first occurs, the 
probability of detection during an encounter can be written as: 
P(N < n) = P(N < m) + P(N = m+1) + ••• + P(N = n) or equally as 

P(N < n) = 1 - [1 - P(N < m)]*(l - (1 - gn) where 

g^ = P(N = i|N < i-1) is the probability of the event detection 

at the i"^^ decision conditioned on the event no detection at an 

earlier decision and 1 < m < n. The second expression is 
generally of greater interest than the first expression, since 
gj^ can usually be more directly related to operational 
parameters such as range and environmental conditions that 
determine a target's detectability than can P(N = i) . 
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with a time rather than the order number to index a decision 
and a random variable T to represent the time index at which 
detection first occurs, P(N < n) becomes P(T < tj^) with 
P(T < tn) = 1 - [1 - P(T < tm)]-[l - g(tn,+i)] ••• [1 - g(tn)] 
where g(tj^) = P(T = tj^|T <tj^_Q^) . 

If g(tj^) << 1 for i = 1,2, ••• ,n, then, to a first 
approximation, ln[l - g(ti)] = -g(tj^) for 1 = 1,2, ••• ,n and 
P(T < tn) = 1 - [1 - P(T < tjjj)]*exp[-2 g(t£)]. This follows 
since P(T < t^) = 1 - [1 - P(T < tj,,) ] •exp[S ln[l - g(ti)] where 
the sum index i = m+1, ••• ,n. A continuous analog to this 

approximation can be used to describe an encounter for which 
g(t£) << 1 for i = m,m+l, ••• ,n and decisions during the 
encounter can be considered to occur continuously. That is, the 
time of an observation corresponding to a decision and the time 
between decisions are both negligible relative to the time of the 
encounter. 

The analog can be developed as follows: First, let St be 

the time between decisions, then tj^ = i*5t and the probability 
of detection P(T < t^) = 1 - [1 - P(T < tm)]-exp[-2: T(ti)*St] 
where t (tj^) = (1/St)*g(ti) is a detection rate function 
(a probability of detection per unit time) and, in terms of St, 
the probability g(t^) = P[T - i*St|T < (i-l)*St]. 

If T is considered to be a continuous random variable, the 
expression for P(T < t^) above indicates that the sum in the 
exponent should be replaced by an integral whose integrand is a 
continuous function T(t). If t ( t) can be determined, then. 
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with g(tj^) as a guide, the cumulative probability of detection 
P(T < t) can be defined by: 



where the limit is for 5t approaching zero. Equation 4 
implies the differential equation: dp(t)/dt = [1 - p(t)]*T(t) 

where p(t) = P(T < t) . A solution to this equation is: 



where t is the time index for a decision during an encounter, tm 
is some time during the encounter and tn > tm. A T(t) that is 
based on a visual detection model is described in Reference 7. 

If the detection capability of a detection system is assumed to 
depend on a target's position relative to the detection system 
during an encounter but not to depend on the clock time, then the 
time index of a decision can be a relative index that determines 
the target position that is associated with a decision rather 
than the clock time associated with the decision. 

The above results apply to the case of an encounter between 
a target and a collection of detection systems. However, if the 
detection systems are not collocated, it is generally convenient 
to describe encounters of this kind in terms of encounters 
between the target and the individual detection systems. In 
either case, if the event target detection for a detection system 
is not independent of the event for other detection systems, then 
in order to describe this in an encounter model the correlation 
between the input to the detection system and the inputs to the 
other detection systems must be specified. This has been done in 



(4) 



T(t) = lim {(l/5t)*P(t < T < t+5tlT < t) ) 



(5) 



'tn 

P(T < tn) = 1 - [1 - P(T < tm) ] 'expC-J tm T(t)dt] 
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some models as follows: First determine the probability of 

detection for each system acting alone. Let be the 

probability that the i^^ system detects the target during the 
encounter under this condition. Next, consider two cases: In 

the first case, the random factors that determine detection for a 
system are independent of those that determine detection for the 
remaining systems. In the second case, the random factors that 
determine detection for the systems are completely dependent. In 
the first case, the probability that at least one system detects 
the target is given by: Pj = 1 - (1 - Pi)*(l - P 2 ) **• (1 ~ Pn) 
where n is the number of detection systems involved. In the 
second case, the probability that none of the systems detect the 
target is given by: 1 - Pq = 1 - Pjjj where Pjjj > Pj^ for 
i = 1,2, ••• ,n since if the m^^ system does not detect the 
target, none of the remaining systems will detect it. The 
probability that at least one system detects the target is given 
by: P = a‘Pp + (1 - a)*Pj where a determines the degree of 
correlation and 0 < a < 1. A way to determine a value for a 
is described in Reference 8. 
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VII. Three Signal Excess Encounter Models 

In the three models described in this section, 
detection is defined in terms of signal excess as it is in 
Section V. Each model determines a cumulative probability of 
detection for a target in an encounter with a passive sonar 
system. An observation in the models is indexed by time and the 
index can usually be considered to be the time at the end of the 
observation. During an encounter, observations are made of a 
series of localization regions. A false alarm can occur for a 
localization region that does not contain a target during an 
observation since the value of RD (or DT) is determined by 
some specified false alarm probability. However, in effect, the 
time to resolve a false alarm is zero. 

To determine signal excess in the models, it is convenient 
to use Equation 2. For each decision in an encounter, there is a 
random variable X(t) defined by Equation 2 that determines the 
random character of the signal excess. For a sequence of 
decisions, the set of these random variables ordered by their 
time index constitutes a stochastic process. And the joint 
distributions of these random variables determines the nature of 
the stochastic process. In the three encounter models described 
in this section, the stochastic process is called a lambda-sigma 
jump process. The time series that are generated by lambda-sigma 
jump processes are represented by the plot in Figure 2 below. 

The jumps in the time series occur at times determined by a 
Poisson process with a mean jump rate lambda. This implies that 
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the time between jumps is a random variable with an exponential 
distribution and that the expected times between jumps t is 
equal to the reciprocal of lambda. The time series that are 
generated by lambda-sigma jump processes are represented by 
Figure 2 below. 




Figure 2. A time series representing a realization of a lambda- 
sigma jump process. On the plot, a in dB equals one unit on 
the vertical axis and t equals one time unit on the horizontal 
axis . 



From Figure 2, note that the observed values of neighboring 
random variables are equal unless a jump has occurred between 
them. When a jump occurs, the first random variable after the 
jump is normally distributed with mean zero and variance a’ and 
it is independent of all the random variables before the jump. 
Conditioned on a jump pattern, this random variable and all the 
random variables between it and the next jump are dependent and 
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the correlation coefficient between any pair is one. That is, if 
the value of the signal excess is known at some time, then all of 
the values between the last jump before that time and the first 
jump after that time are also known. However, since the jumps 
occur randomly, knowing the value of the signal excess with 
certainty at some time does not determine the values of the 
signal excess with certainty at neighboring times. In the 
unconditioned case, the correlation coefficient between the 
random variables X(t) and X(t+T) is equal to 1/e. For this 
reason, t is referred to as a relaxation time. 

It appears that the use of the lambda-sigma jump process is 
based more on past practice than on experimental justification. 

In this regard, see Reference 9. By referring to Equation 1, it 
can be seen that the lambda-sigma jump process is determined by 
the sum of the stochastic processes that determine the random 
variables on the right side of this equation. Although the sum 
of a collection of normal random variables is a normal random 
variable, the sum of a collection of lambda-sigma jump processes 
is not a lambda-sigma jump process. This suggests that if the 
lambda-sigma jump process does adequately describe the 
variability of the signal excess, then the majority of the 
variability of the signal excess may be due to a single one of 
its elements. For example, transmission loss. 

In the three encounter models described below, detection is 
defined in terms of signal excess as described in Section IV and 
decisions are indexed by a time that can usually be considered to 
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be the time of the decision. During an encounter, observations 
are made of a series of localization regions. For localization 
regions that do not contain a target, the signal observed during 
an observation of the region is zero. In the models, the time to 
resolve a false alarm is zero. However, false alarms are not 
ignored in that the value of RD (or DT) is determined by some 
specified false alarm probability. Consequently, although the 
time to resolve a false alarm is zero, the cost associated with a 
false alarm is not zero. 

The First Passive Sonar Encounter Detection Model: This 

model describes an encounter in terms of a series of decisions 
with each decision based on the signal excess X 3 j;(t) at a time 
corresponding to the end of an observation. The observations are 
of equal duration and the integration time that determines the 
recognition differential is equal to the duration of the 
observations. In the model, X 3 g(t) is determined by a lambda- 
sigma jump process. For an encounter in which SE(t) is 
unimodal and in which the time of the single maximum is prior to 
the end of the encounter, it is shown in Reference 10 that the 
probability p that detection will occur during the encounter 
is given by the following equation: 

(6) p = 1 - [(1 - Pc)/(1 - P-Pc)]*(l - P-Pi) ••• (1 - P*Pm) 

where p = 1 - exp(-St/T) and p^ = $[SE(tj^)/a] for i = 1,2, 
••• ,m. Here, 5t indicates the duration of an observation and 
i indicates the standard normal cumulative distribution function 
as before. The integer c is the index of the decision time t^ 
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for which SE(t<;;) is greater than or equal to SE(tj^) for any 
time tj^ and ti < tc < tin- 

As T approaches zero, p approaches one and Equation 6 
approaches this form: 

(7) p = 1 - (1 - Pi) ••• (1 - Pm) . 

In this limit, the signal excess random variables are all 
independent. Note that Equation 7 applies without the condition 
that SE(t) be unimodal. 

As T approaches infinity, p approaches zero and 
Equation 6 approaches this form: 

(8) p = Pc- 

In this limit, the correlation coefficient between any pair of 
signal excess random variables is equal to one. Note that 
Equation 8 applies without the condition that SE(t) be 
unimodal. Equation 8 defines a complete dependence encounter 
model . 

The Second Passive Sonar Encounter Detection Model: This 

model is in a sense a third limiting form of the first passive 
sonar encounter detection model. In this limit, the time between 
decisions approaches zero. However, in this limit the 
integration time that determines the recognition differential is 
not equal to St and it does not approach zero. It is, in 
effect, chosen by the user of the model through the user's choice 
of the value for the recognition differential. For an encounter 
that begins at tl and ends at t2 and for which XsE(t) is 
determined by a lambda-sigma jump process and SE(t) is 
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unimodal, it is shown in Reference 10 that for this limit, 
Equation 6 has the following form: 



■t2 

(9) p = 1 - [1 - p(tc) ] •exp[-(l/T) .Jti p(t) dt] 
where p(t) = #[SE(t)/a] and where now t^, is the encounter 
time such that SE(tc) is greater than equal to SE(t) for any 
other encounter time t and tl < t^ < t2 . 

The Third Passive Sonar Encounter Detection Model: This 

model describes an encounter between a target and a passive sonar 
detection system in which detection occurs during an encounter if 
the average value of the square of the continuously observed 
signal-to-noise ratio over a time interval of length u is 
greater than or equal to the square of the signal-to-noise ratio 
that determines the recognition differential for an integration 
time equal to u. With R(s) the random signal-to-noise ratio 
at a time s and random signal-to-noise ratio that 

determines the random recognition differential for an integration 
time u, detection during an encounter occurs at the first time 
t that the following inequality is satisfied: 



(10) (1/u) Jt-u [R(s)/Rm(u) ] * ds > 1 

where the time origin is chosen so that t > 0 and where the 
integration time u = t for t < to and u = to for t > to 
where to is a maximum integration time. The random integrand 
in the inequality is related to the random signal excess at the 
time s for an integration time u. The relation is: 

(11) 10 log [R(s)/R„,(u) ] * = 2[SE(s;u) + X(s)] 
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where SE(s;u) is the expected value of the signal excess at a 
time s for an integration time u and X(s) is the random 

component of the signal excess at the time s. In the model, 

X(s) is determined by a lambda-sigma jump process and SE(s;u) 
is determined by an expected value sonar equation with a 
recognition differential RD(u) = 10 log rju(u) . Here, Tjjj(t) is 
the value of the signal-to-noise ratio that gives a probability 
of detection equal to . 5 for an integration time t and a 
specified probability of false alarm Pf. With the signal 
detection process described by a Case II signal detection model, 
the detection index necessary to give the required operating 
point (pf,.5) is related to the integration time t and the 
signal-to-noise ratio rm(t) by: 

(12) d = u- (BW) [ri„(t)]* 

where BW is the bandwidth of the receiver. For a spectrum 
analyzer, BW would be the bandwidth corresponding to a given 
frequency resolution and d would be the detection index 
required in order to be at the operating point (Pf/.5) for a 
signal that was contained within a bandwidth BW. Since d in 
Equation 12 must be the same for t = u and t = tQ, 

(13) RD(u) = 5 log(to/u) + RD(to) 

where to is the maximum integration time. Then, since 
SE(s;u) - SE(s;to) = RD(to) - RD(u) , by using Equation 13 and 
Equation 11, Relation 10 becomes: 

ft (1/5) [X(s) + SE(s;to) - 5 log(to)] 

(14) Jt-u 10 ds > 1 
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where as above the time origin is chosen so that t > 0, the 
integration time u = t for t < to and u = to for t > to 
and where SE(s;to) is the expected value of the signal excess 
at the time s for a recognition differential determined by an 
integration time to* In an encounter, detection occurs the 
first time that Relation 14 is satisfied. 

As is pointed out in Reference 11, the appeal of the Third 
Passive Sonar Encounter Detection Model relative to the Second 
and First Passive Sonar Encounter Detection Models is that it 
appears to more closely describe the detection process in passive 
sonar detection systems that display their processed data to an 
operator in a continuous manner over a time window of duration 
tQ. However, results reported in Reference 12 indicate that the 
difference between the three models may not be significant in 
some types of encounters. 
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VIII. straight Line Encounters 

Suppose a target's detectability depends on its range 
from a detection system and that the probability of detection is 
effectively zero beyond a range rj„ for any target azimuth. In 
this report, an encounter between the target and the detection 
system is the event that the range between the target and the 
detection system is less than or equal to r^. In addition, 
suppose rjj, is small enough so that when the target and the 
detection system are having an encounter they can be considered 
to be moving on planes parallel to a tangent plane to the earth's 
surface at some point in their vicinity. If this is the case, 
then while the target and detection system maintain a constant 
course and speed during an encounter, the encounter is called a 
straight line encounter. 

A straight line encounter can be described in terms of a two 
dimensional rectangular coordinate system whose plane is parallel 
to the tangent plane to the earth. If the coordinate system is 
stationary relative to the detection system with the detection 
system located at the origin and is oriented so that the target's 
motion is parallel to the y-axis and is in the positive 
y-direction, then the target's x-coordinate during a straight 
line encounter will be constant. The constant is equal to the 
target's horizontal range at the closest point of approach (CPA) 
on the straight line track on which the target is moving relative 
to the detection system during the encounter. This range is 
called the target's lateral range. 
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A complete straight line encounter is a straight line 
encounter that begins at a range from a detection system that is 
greater than or equal to rjj^ and continues past CPA to a range 
from the detection system that is again equal to or greater than 
r^. Let p(x) be the cumulative probability that a target is 
detected by a detection system in a complete straight line 
encounter in which the target's lateral range is x. Then the 
function p(x) defines what is called a lateral range cuirve or 
lateral range function. 

Let p be the probability that a target is detected during a 
complete straight line encounter. If the lateral range of a 
target in a straight line encounter is assumed to be a continuous 
random variable X with a uniform distribution with fx(^) “ 
for |x| < a/2 and p(x) = 0 for |x| > a/2, then the 
probability that a target will be detected during a complete 
straight line encounter is given by: 



where the limits of integration can be used since the value of 
p(x) is zero for |x| > a/2. Equation 15 suggests a measure of 
a detection system's capability to detect a target in a straight 
line encounter. The measure W is called sweep width and 



(15) 




dx 



(16) 




dx. 



With this definition. Equation 16 becomes: p = (l/2a) *W. 
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IX. Two Intermittent Signal Encounter Models 

In the intermittent signal encounter models that are 
described in this section, an encounter is a complete straight 
line encounter, and during an encounter a target either emits a 
signal (an acoustic transient) or its presence (a visible 
submarine mast) is the cause of a signal at various times. Two 
cases are considered: In the first case, the signals occur 

periodically, the signals are of length 6t and the time between 
the occurrence of signals is t where t > St. In the second 
case, St = 0 (the signals are instantaneous) and the signals 
occur at times determined by a Poisson process for which the 
expected time between signals is equal to t. In the model, the 
detectability of a target signal depends on a target's horizontal 
range from a detection system, but on no other factors. If a 
signal occurs while the target is within a range r, it will be 
detected. For a continuous signal, the lateral range function of 
a detection system for a target is: p(x) =1 for |x| < r and 
p(x) = 0 for |x| > r where the horizontal range r is 
determined by the characteristics of the detection system and the 
target. The geometry for an encounter is shown in Figure 3 
below. 

For intermittent signals, the length of a target's track 
relative to a detection system on which a signal will be detected 
is 2»(r* - x*)5 + wSt where w is the speed of the target 
relative to the detection system. So, a target's exposure time 
during an encounter is (2/w) • (r* - x*)^ + St. 
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For periodic signals, there are two cases. In the first 
case, r > w* (t - 6t)/2. In this case, the signals result in the 
following lateral range function: 
p(x) = 0 for |x| > r 

(17) p(x) = 1 for |x| < {r^ - [w- (t - fit)/2]M^ 
p(x) = [2/ (w T) ] • (r^ - x^ ) ^ + St/T otherwise 

In the second case, r < w(r - 5t)/2 and the middle equality 
in Equation 17 does not apply. 

For signals that are instantaneous and whose occurrence is 
determined by a Poisson process, the signals result in the 
following lateral range function: 

p(x) = 1 - exp { - [ 2/ (w* T ) ] ♦ (r^ - x^ ) ^ } for |x| < r 

(18) 

p(x) = 0 for |x| > r. 




Figure 3 . The encounter geometry for the two intermittent signal 
models described here. 
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For signals whose occurrence is determined by a Poisson 
process and for which St > 0, signals can overlap. If this is 
allowed, then Equation 18 can be modified to describe this case 
by adding St/T to the term in the exponent of Equation 18 that 
is within the square brackets. In particular, note that this 
modified Equation 18 can be approximated by the bottom equality 
in Equation 17 when (2/wT)(r* - x*)= + St/T << 1. This 
implies that when the expected time t between signals is large 
relative to the exposure time (2/w) (r® - x®)^ + St, the periodic 
occurrence model and the random occurrence model are effectively 
equivalent. 
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X. A Random Search Model 



A search of a region in which a target moves on a track 
that consists of a number of straight segments placed in such a 
way that in a limiting sense every section of the region is 
equally likely to be searched on a segment is referred to as a 
random search in Reference 6. Representation of a search region 
with the track segments that could be imagined to be the tracks 
of a random search are shown in Figure 4 . 




Figure 4 . A search region and a track that could be described as 
a random search track. 

Two developments of a model to describe this kind of search 
are contained in this section. The first development is based on 
the following conditions: 1. A target is at a fixed position 

within a defined search region. 2. A searcher's track is a 
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sequence of straight line segments that are within the search 
region. 3. The searcher's detection system is such that while on 
a track segment, a rectangle is searched that is contained within 
the search region, is of length equal to the length of the track 
segment and is oriented so that its long axis is parallel to the 
track segment. 4. The probability that the searcher's detection 
system will detect a target while on a track segment with a 
search rectangle that does not contain the target is zero. The 
probability that the searcher's detection system will detect a 
target while on a track segment with a search rectangle that 
contains the target is p(x) where x is the target's lateral 
range for the track segment and p(x) is the lateral range curve 
for a complete straight line encounter lateral range x. A 
representation of a search rectangle is shown in Figure 5 below. 
5. The track segments are located in such a way that the event 
that the target is within the search rectangle associated with a 
track segment is independent of the event that the target in the 
search rectangle associated with any other track segment. And 
the probability of the event is equal to the ratio of the area of 
the search rectangle to the area of the search region and, given 
a target is within a search rectangle, its position is uniformly 
distributed over the rectangle. 

Condition 4 implies that the random search model is based 
on the concept of a complete straight line encounter. The 
definition of an encounter that is intended here is that given in 
Section VI. This implies that in the random search model the 
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time to resolve a false alarm is zero. However, for the model, 
and pf are considered to be determined by some criterion 
such that Pf is less than one. Consequently, although the time 
to resolve a false alarm is zero in the model, the cost 
associated with a false alarm is not zero. (A concise model that 
accounts for the time to resolve false alarms is described in 
Reference 13.) Condition 4 also implies that when a searcher is 
on a track segment with a search rectangle that contains a 
target, the encounter is a complete straight line encounter. And 
Condition 5, which can be considered to specify a random 
arrangement of the track segments, implies that when this is the 
case, for the complete straight line encounter, the target's 
lateral range is a random variable that is uniformly distributed 
between -b/2 and b/2 where b is width of the search 
rectangle (the dimension of the rectangle perpendicular to the 
associated track segment) . 



Swept Area 



Track 

Figure 5. A track segment and its associated search rectangle, 
that could correspond to a search with an aircraft mounted 
infrared detection system. 
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Based on the above considerations, the probability that a 
target will be detected while a searcher is on a track segment 
with an associated search rectangle that contains the target is 
given by: 



where fx(^) “ -b/2 < x < b/2 and fx(^) “ ® 

p(x) = 0 otherwise. Note that the left side of Equation 19 
applies to any complete straight line encounter in which the 
target's lateral range for the encounter is considered to be a 
random variable with a distribution determined by the probability 
density function fx(^) ♦ given that the target is 

within the search rectangle associated with a track segment, then 
the unconditional probability that the target will be detected on 
the track segment is given by: (W/b) • (SA/A) where SA is the 

area of the search rectangle associated with the track segment 
and A is the area of the search region. With 1 the length of 
the rectangle, 5 A = b*l and the probability becomes: (W*l)/A. 

Then, since the event that the target will be in the search 
rectangle of a track segment is independent of the event that it 
will be in the search rectangle of any other track segment, the 
probability p that a random search consisting of m track 
segments will detect the target is given by: 

1 - [1 - (W*li)/A][l - (W-l 2 )/A] ••• [1 - (W*ln)/A] where 1^ 
is the length of the i^^ track segment. The probability is 
also given by: p = 1 - exp{E ln[l - (W*lj^)/A]} where the sum 



(19) 




fx(x) dx = W/b 
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index i = 1,2, ••• ,n. If (W*li)/A << 1 for i = 1,2, ••• ,n, 

then this expression can be approximated by: 

(20) p = 1 - exp[-(W*l)/A] 

where 1 = S is the track length of the search. Equation 20 

is known as the random search formula. 

The second development of the random search formula is based 
on Equation 5 and a detection rate for a random search given by: 

T (t) = W*v(t)/A. With this detection rate and Equation 5, the 
random search formula is given by: 

(20a) P(T < t) = 1 - exp ( - [W- 1 (t) ]/A) 

where l(t) is the track length for a random search that starts 
at time 0 and ends at time t and 



't 

(20b) l(t) = ^0 v(s)ds. 

Replacing P(T < t) by p and 1 (t) by 1 gives Equation 20. 
In the form of Equation 20a, the random search formula indicates 
explicitly the relation between the probability of detection and 
the duration of a random search. Note that Equation 20a implies 
that the sweep width is independent of speed over the range of 
speeds in the encounter. 

Reference 14 contains an example of an application of the 
method used in the second development of the random search 
formula to a random search where the search region expands with 
time. 
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XI. Ladder and Barrier Search Models 



In some barrier searches, the barrier search track is a 
ladder search track relative to a reference system that moves 
with the target. This fact is used in the barrier search model 
development that follows the two ladder search model developments 
below. The first ladder search model is referred to as an ideal 
ladder search model because of the idealizations that are 
involved in its description of a ladder search. The second 
ladder search model is referred to as a degraded ladder search. 

It can be considered to describe a ladder search track in which 
navigational errors result in omissions and overlaps in coverage. 

An Ideal Ladder Search Model: The model is based on the 

following conditions: 1. A ladder search region is a rectangle 

that contains a fixed target. 2. During a search of the region, 
the searcher's detection system moves on a set of m parallel 
track segments of length b separated by a distance s. 3 . As 
the detection system moves along a track segment, it searches a 
rectangular strip of length b and width s within the search 
region. 4. The m rectangular strips that correspond to the m 
track segments completely cover the ladder search region with no 
overlap. 5. If a target is within the rectangular strip 
corresponding to a track segment, then there will be a complete 
straight line encounter between the target and the detection 
system when the detection system moves along the track segment 
and the lateral range of the encounter will be uniformly 
distributed across the width of the strip. If the target is not 
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in the rectangular strip, then there will not be an encounter and 
the probability that the target will be detected while the 
detection system is on the track segment is zero. 




Figure 6. A schematic representation of a ladder search geometry 
for a case in which the ladder search track segments are 
superimposed on and bisect their corresponding rectangular 
strips. 

Since targets outside of the rectangular strip that 
corresponds to a track segment cannot be detected while a 
detection system is on the track segment because of Condition 5, 
in the model, the sweep width W of a searcher's detection 
system must satisfy the relation W < s . In particular, W = s 
only holds when the detection system detects a target that is in 
a rectangular strip with probability one for any target lateral 
range. This kind of detection system is sometimes referred to as 
a cookie cutter detection system. However, this terminology can 
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be misleading since it suggests the detection system detects 
equally well for all azimuths. But this is not a requirement on 
the system in order that W = s. 

The ideal ladder search model implies that if the conditions 
of the model are satisfied, then the probability p that a 
target will be detected by a an ideal ladder search is given by: 

(21) p = W/s 

where W/s < 1. The quantity W/s is called the coverage 
factor in this case. 

A Degraded Ladder Search Model: The above model implies 

perfect navigation in addition to other idealizations. A model 
of a ladder search is given in Reference 6 that could be used 
for cases in which this is a poor assumption. The model which is 
referred to here as a degraded ladder search model can be 
considered to describe navigational inaccuracies in terms of 
omissions and overlaps of the rectangular strips. It can be j 

developed as follows: Consider a random search in the ladder 

search region whose track length is equal to the search track 
length required to complete an ideal ladder search, that is, a 
track length 1 = m*b. The degraded ladder search model 
describes the result of omissions and overlaps in a ladder search 
to be such that the probability of detection for this random | 

search is equal to the probability of detection for the degraded , 

ladder search. Consequently, since the area of the ladder search j 

I 

region is m*s*b, for the degraded ladder search model: 

i 

(22) p = 1 - exp (-W/s) . l’ 
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Here, the requirement that the coverage factor W/s < 1 for 
Equation 21 can be relaxed. However, it should still be 
considered as an approximate condition. 

The condition that the target be fixed within the 
rectangular search region is critical to both Equation 21 and 
Equation 22. However, these results are also applicable to a 
search for a moving target under the conditions that are 
described next. 

A Barrier Search Model: A target moves with a constant 

course and a constant speed u. Both the target's course and the 
target's speed are known by a searcher. The searcher establishes 
a barrier of width b that is perpendicular to the target's 
track and moves on the barrier with a speed v > u. The barrier 
is designed so that in a reference system relative to the target 
the barrier search is a ladder search that satisfies the 
conditions for a ladder search that are given above. There are 
two cases to consider: 1. The barrier is established in front of 
the target. 2. The barrier is established behind the target. 

From the search geometry for a barrier established in front 
of the target, it can be seen from Figure 7 below that 
0 = sin~^(u/v) and d = vt where t = s/(v + u) is the time to 

move from one search leg to the next. The angle 0 and the 

perpendicular distance d that depend on u, v and s, and the 

width of the barrier b are the quantities that are required in 
order to establish the barrier operationally. 
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b 

Figure 7. A barrier search track shown for a barrier established 
in front of the target. The track is shown in a reference systein 
fixed relative to the earth. 

For a barrier that is established in front of a target, one 
of three barrier types will result. A barrier's type is 
determined by the relation of the distance d to the distance 

g = ut where the time t = b/ (v^ - u^ ) ^ is the time to complete 

a search leg (cross the barrier) . The barrier type is determined 
as follows: 1. For g < d, the barrier is an advancing barrier. 

2. For g = d, the barrier is a stationary barrier. 

3. For g > d, the barrier is a retreating barrier. 

For a barrier established behind the target, there is only 
one barrier type and it is called an overtaking barrier. For an 
overtaking barrier, 0 = sin“^(u/v) as for a barrier established 
in front of the target. But, for an overtaking barrier, 
r = s/(v - u) and d = vs/(v - u) . 
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Given the target crosses the barrier, the probability of 
detection for an ideal barrier search is given by Equation 21 
and the probability for a degraded barrier search is given by 
Equation 22 where the terminology refers to the nature of the 
ladder search in the reference system moving with the target. 

A discussion of an application of these two equations to a search 
for a magnetic anomaly target is given in Reference 15. 
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XII. A Target State Estimation Procedure 

A target state estimation procedure based on bearing 
observations is developed in this section that generates point 
estimates of a target's position and velocity vector coordinates 
in a rectangular coordinate system. The procedure is based on a 
model in which bearing errors are unknown and are not determined 
by random variables with known distributions. Because of this, 
confidence regions for the estimates are not generated by the 
procedure. However, for a moving target, it illustrates general 
characteristics of bearings only target motion analysis (TMA) . 

The model is defined as follows: 1. The target moves in a plane 

with a constant but unknown course and speed. 2. Observations of 
the target are made from known positions at known times. 3_. The 
observations provide only target bearings with unknown errors. 

The model geometry is shown in Figure 8. 




Figure 8. The geometry of the target motion analysis model. 
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The procedure criterion is: For observations from n positions, 

choose target position estimates and target velocity component 
estimates Ujj and Uy that make the sum of the squares of the 
algebraic distance between the estimated positions and their 
corresponding observed bearing lines a minimum. From Figure 8, 
it can be seen that the algebraic distance can be written as 
dj^ = [x-t(i) - Xo(i)]*cos 0^ - [Yt(i) ~ yo(i)]*sin 0j^. Because of 
the requirement that the target move with constant course and 
speed during the encounter, the number of independent estimates 
is reduced from 2n to 4, Ujj / Uy and any two position 
estimates x-t(j), yt(j)* "the following development, j = 1 

and with i = 2,3, ••• ,n the remaining estimates are given by: 

xt(i) = xt(l) + Ux*(ti - ti) and yt(i) = Yt(l) + ^y* (tf - t^) . 
To determine "best” estimates of the target state parameters, 
take the partial derivative of the sum S = S (dj^) * with respect 
to each of them. Then set the four partial derivatives equal to 
zero. This creates four linear equations in x-t(l) / Yt(i)f 
and Uy whose solution are the desired estimates x^{l) , Yt(^)' 
Ujj and Uy. In matrix notation, the equations can be 
represented by AX = B where the elements of X are: 
xii = xt(l) , X 21 = yt(l)/ ^31 = '^x ^41 = '^y* ^ necessary 

condition for a unique solution for X is that n > 4. 

Otherwise, the determinant of A will be equal to zero. The 
procedure can also be used if a target's course and speed are 
constant and known and, in particular, if the target is 
stationary so that Uj^ and Uy both equal zero. In this case. 
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since the number of unknowns is two, the number of linear 
equations is also two and a necessary condition for a unique 
solution is n > 2. 

Now, suppose the observations are at positions and times 
that correspond to the positions and times of an observer moving 
on some constant course at some constant speed (including zero 
speed) . In this case, the observation position coordinates are 
related by the following equations: XQ(i) = Xq(1) + ~ '^l) 

and yo(i) = Yo(^) ''^y('^i “ where Vj^ and Vy are the 

required velocity components of the observer. Using these 
equations of motion, the matrix equation AX = B can be 
transformed to the matrix equation Ax'= 0 where the elements 
of the matrix X are related to the elements of the matrix x' 
by the equations: x^i = Xt(l) - Xq(1), X 21 = yt(l) - YoCl)' 

^31 “ ^x " "'^x ^41 “ ^y ” 

Since the linear equations represented by Ax' = 0 are 

homogenous, they do not have unique solutions and consequently 
neither do the equations represented by AX = B. However, if 
there is at least one observation whose time and position is not 
determined by the above equations of motion, then the 
transformation from X to x' cannot be made, and in general a 
unique solution for X can be found. If the observations are 
made from a platform that is moving with a constant course and 
speed, this condition can be achieved by either changing the 
course, the speed or both prior to completing the observations. 

Estimation models that describe bearing error as a random 
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variable provide a basis for determining confidence regions for 
point estimates. A model is developed in Reference 16 that does 
this for either target bearing observations made from two or more 
points simultaneously or for a target that is stationary relative 
to the observation points. 
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XIII . Position Distributions That Change with Motion 

Target motion models provide a basis for determining 
position distributions that change with target motion. In this 
section, two classes of target motion models are considered. In 
the first class, a target moves in a plane with a constant course 

and speed. In the second class, a target moves in a plane but 

its course or speed changes during the motion. Three members of 
the first class are developed first. This is followed by a brief 
discussion of some models of the second class. 

The First Motion Model: In the first model, X(0) and 

Y(0), a target's random position coordinates at time zero, have a 

joint distribution that is circular normal with mean vector 

(0,0) and standard deviation a. The joint distribution of the 
velocity components Ujj and Uy is circular normal with mean 
vector (Ujj,Uy) and standard deviation and the random 

variables Ujj and Uy are independent of X(0) and Y(0). The 
target coordinates at time t are given by X(t) = X(0) + Uj^*t 
and Y(t) = Y(0) + Uy*t. Therefore, the random vector 
[X(t),Y(t)] is a normal random vector and its characteristic 
function is equal to the product of the characteristic functions 
of the random vectors [X(0),Y(0)] and [Uj^‘t,Uy‘t] . This 
implies that X(t) and Y(t) have a circular normal 
distribution whose mean vector is (Uj^*t,Uy*t) and whose 
variance is (a* + a^*t*). For a more complete discussion of the 
basis for the above argument, see Reference 7 
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The Second Motion Model: In the second model as in the 

first, X(0) and Y(0), a target's random position coordinates at 
time zero, have a joint distribution that is circular normal with 
mean vector (0,0) and standard deviation a. However, in this 
model, target speed u is known, but its course is equally 
likely to have any value between 0 and 27 t, so Ujf = u*sin # 
and Uy = u*cos where $ is a random variable with a uniform 

distribution and 0 < 0 < 2tt. In this model, only the random 
variable § is required to determine the target's velocity. The 
equation that is described below that determines the distribution 
of the random position coordinates [X(t),Y(t)] can also be used 
to determine it for the first model. It is based on the 
following consideration: If a target's course and speed do not 

change during its motion, the joint density function of its 
position coordinates at some time t is determined by the 
equation: 



(23) fx(t) ,Y(t) (^'Y't) = oofx(O) , Y(0) (q'S;0) fv,w(v,w) dvdw 

where q = x - vt and s = y - w*t and V = Uj^, W = Uy, v = Uj^ 
w = Uy and t > 0. Equation 23 can be developed as follows: 

To first order, fx(t) , Y(t) probability that 

a target's coordinates are in an element of area 5x*Sy and for 
given values of v, w and t, target positions in an element of 
area Sq-Ss will be translated to an element of area Sx*6y 
that is identical in form and size to 5q*Ss. And so, to first 
order, fx (o) / y ( 0 ) (< 3 ' ^ * ^V, W v) • Sv Sw is the 

probability that the target's coordinates at time 0 are in an 
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element of area 5q*5s that is located such that the target's 
coordinates will be in the element of area 5x*6y at time t 
since x = q + vt and y = s + w*t. And, to first order, the 
sum of such probabilities for all pairs of values of v and w 
is also the probability that the target's coordinates at time t 
are in the element of area 6x*6y. In the limit after equating 
the two expressions for this probability and cancelling the 
common factor 5x*5y, Equation 23 results. 

With the velocity components in polar coordinates, the 
integral of Equation 23 is a single integral over <p and the 
integrand of the integral is (l/2Tra*) •exp[-(q* + s*)/2a*] • (1/2 tt) 
where now q = x - u*t*sin 0 and s = y - u*t*cos 0. The result 
of the integration is: 

(24) (l/2Tta*) exp{-[x* + y* + (u*t)*]/2a*) Io[(x* + y*)5.u*t/a*] 

where t > 0 and Iq indicates the hyperbolic Bessel function of 
zeroth order. In Reference 6 for several values of t. 
Expression 24 is plotted in terms of r = (x* + y*)^, the 
target's range from the origin. The plot shows a limiting 
characteristic of this distribution that can be indicated as 
follows: First replace (x* + y*)^ by r in Expression 24. 

Then, by first multiplying and then dividing the resulting 
expression by exp (-r *u* t/a* ) ' the following expression results: 

(25) 1/ (2Tra* ) exp{ - [ 1/ (2a * ) ] (r - u* t) * } Iq (r*u* t/a * ) exp (-r *u • t/a* ) 
where t > 0. The second factor in Expression 25 is 
proportional to the density function of a normal random variable 
whose mean is u*t and whose standard deviation is a. And the 
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last two factors are such that their product decreases slowly as 
a function of r. The consequence of this is that a plot of 
Expression 25 against r for values of t greater than 4*a/u 
has the appearance of a normal density function. In particular, 
this indicates that u*t in the exponent of the second factor in 
Expression 25 can be viewed as the radius of an average furthest 
on circle. 

A target's random rectangular coordinates X(t) and Y (t) 
and its random bearing 0(t) and range R(t) from the origin 
are related by: X(t) = R(t) sin 0(t) and Y(t) = R(t) cos 0(t). 

By using these equations, Expression 24 can be transformed to: 
(26) ( 1/2 it) (r/2TTo ’ ) exp{ - [r * + (ut) ’ ]/2a’ }Iq (rut/o’ ) 

which is the joint density function of the random variables R(t) 
and 0(t). The marginal density function of R(t) can be 
obtained from Expression 26 by integrating this joint density 
function over the possible values of 0(t) which in radians is 
over the interval 0 to 2 tt. Consequently, the marginal density 
for R can be obtained from Expression 26 by multiplying it by 
2 tt. Tabulated values of the cumulative distribution function 
Fj^(r;t) for the marginal distribution of R for the second case 
are listed in Reference 17. 

The Third Motion Model: In the third model, only a target's 

maximum speed Uj„ is known, but a target's position at time zero 
is known to be at the origin of the position coordinate system 
and, after that time, to be uniformly distributed on a circular 
disk of radius Ujj,*t centered at the origin. This implies that 
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(27) 


fx(t) ,Y(t) (X/Y'-t) 


= l/(u^-tM 


for 


t > 0 where x* + y* 


< u^*t* and that 


(28) 


^R(t) ,0(t) (r-fCt) = 


^ r/(7r-u^*t*) 


for 


t > 0 where 0 < r < 


Uju*t and 0 < a 



The range of values of r and a are independent and the joint 
density function of R(t) and 0 (t) is equal to the product of 
1/(2*7 t) and 2*r/(u^*t*). This implies that for the marginal 
distributions with density functions: iR(t) “ 2 -r/(u^.t*) 

where 0 < r < Ujjj*t and f 0 (-t) “ l/( 2 *ir) where 0 < a < 2 *tt, 

the random variables R(t) and 6 (t) are independent, their 
joint density function is given by Equation 28 and the joint 
density function of X(t) and Y(t) is given by Equation 27. 
These two marginal distributions can be achieved by choosing, at 
time zero, a course © from the uniform distribution with 
density function f 0 (a) = 1/(2 •it) where 0 < a < 2 *‘ 7 r and a speed 
U from the triangular distribution with density function fuC^^) 
= 2 *u/u^ where 0 < u < Uj„. 

Motion Models of the Second Class: In the second class of 

motion models, a target's course or speed or both can change. In 
general, a monte carlo simulation method is required in order to 
determine a position distribution that is based on such a model. 
As an example of cases in which the distribution can be described 
analytically, see Reference 18 and Reference 19. As an example 
of a case in which it can not, suppose a target's initial 
position is described in terms of a number assigned to a 
subregion in the xy-plane where the number assigned represents 
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the probability that the subregion contains the target at an 
initial time. In addition, suppose for each subregion a course 
and speed distribution is determined by assigning numbers to 
course and speed pairs where a number represents the probability 
the target will have the course and speed at the initial time 
given it is in the subregion at that time. Next suppose for each 
course and speed pair there is a time distribution that 
determines the duration of the course and speed pair and that the 
time distribution is determined by a number assigned to each 
discrete time point where the number represents the probability 
that the target's course and speed pair will be determined by a 
new course and speed distribution. By extending this kind of 
procedure and then implementing it in a monte carlo simulation, 
one can generate complex position distributions that describe a 
target's position at discrete time points. 
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XIV. Position Distributions That Change with Search 

Suppose a target's position at some time is described 
by a position distribution. Now suppose information becomes 
available that a search has been conducted for the target and 
that the target has not been detected. Or suppose the 
information is that the target has been detected. In the first 
case, negative information is available that can be used to 
modify the position distribution. In the second case, positive 
information is available that can be used to modify the position 
distribution. In both of the cases, a target's position 
distribution is assumed to be specified in terms of a set of 
discrete probabilities where each probability corresponds to a 
subregion of the region that contains the target and each is the 
probability that the target is in that subregion. 

Position Distributions and Negative Information: For a 

region that contains a target and consists of n subregions, let 
the event = {the target is in the i^^ subregion}. And let 

the event C = {no contact}. Then, given no contact in a search 
of the region, the targets 's position distribution can be 
modified as follows: 

(29) P(Si|C) = P(c|Si) •P(Si)/P(C) 

where i = 1,2, ••• ,n and P(C) = 2 P(C| Sj ) • P(Sj ) with the sum 

index j=l,2, ••• ,n. Note that Equation 29 can be obtained 

by using Bayes theorem. To illustrate how Equation 26 might be 
used, suppose that a search in a subregion is considered to be a 
random search and that the sweep width of detection system 
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against the target depends on subregion being searched. For this 
case, let be the area of the i^^ subregion and let be 

the sweep width in that subregion. Then, given no contact in a 
search of a subregion, P(c|Sj^) = exp[-(Wj^*l£)/Aj^] where Ij^ is 
the track length of the searcher in the i^^ subregion. Given 
values for P(Sj^), Vl^/A^ and Ij^ for i = 1,2, ••• ,n, a 
position distribution can be determined that has been modified by 
the negative information. 

Position Distributions and Positive Information: In the 

case of positive information, the event C = {a contact) occurs. 
This event is the union of two mutually exclusive and exhaustive 
events: T^ = {a true contact) and = {a false contact). 

Relative to the Venn diagram of Figure 1, the event defined 

above corresponds to C corresponds to T^, 

corresponds to (D^ n H^) and F^ corresponds to (D^ n Hq) . 

For generality, suppose true contacts do not localize a target to 
a single subregion. Then after a search of a region that has 
resulted in a contact, the target's position distribution can be 
modified as follows: 

(30) P(Si|C) = P(Si|Tc)*P(Tc|C) + P(Si|Fc) -PCFcIc) 

where i = 1,2, ••• ,n since T^ = T^ o C, F<^ = F^ n C and 

P(Si|c) = {P[(S£ n Tc) + P(Si n Fc)])/P(C). The probability 

p = P(TclC) has been called the credibility of the contact. In 
terms of p Equation 30 becomes: 

(31) P(Si|C) = P(Si|Tc)-p + P(Si|Fc)-(l - P) . 
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For the model, PCSj^ITq) can be determined by the coverage 
characteristics of the detection system that is used to make the 
contact. In particular, PCS^ITq) = P(Tcl S^) •P(Si)/P(Tc) where 
P(Tq) = S P(Tq I S j ) • P(Sj ) and the sum index j = 1,2, ••• ,n. 

Then, using the correspondence between the events and 

and the events Tq and (D^ n , a definition for P(Tc|Si) 
is P(TQ|Sj^) = (Pd) i where (Pd) i average probability of 

detection over the i^^ subregion. Note, if C = T^, then p = 1 
and Equation 31 can be written as P(Si|C) = P(C|Si)/P(C) where 
i = 1,2, ••• ,n and P(C) = S P(C | S j ) • P (S j ) with the sum index 
j = 1,2, ••• ,n which is analogous to Equation 30. In one 
positive information model, P(Si|F^,) = P(Si) (false contacts 
supply no information about a target's location), and p is 
determined subjectively based on factors associated with the 
detection system used to make the contact. The probability 
determined by 1 - p has been called the false alarm 
probability. However, it is not the probability pf = P(D 2 ^|Hq). 
Since 1 - p = P(F<^|C), it is the probability P(Hq|Di). 

This approach could also be used to modify a position 
distribution given negative information. To do this, note that 
the event C is the union of two mutually exclusive and 
exhaustive events: the event T^ = (a true no contact) and the 
event F^ = {a false no contact). Relative to the Venn diagram of 
Figure 1, C corresponds to Dq, Tq corresponds to (Dq n Hq) 
and Fq corresponds to (Dq n and the countrpart of 

Equation 31 is: 
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(32) P(Si|C) = P(Si|Tc)-p + P(Si|Fc)-(l - P) 

where p = PCT^Ic). Relative to the Venn diagram of Figure 1, 
p = P(Do|Hq). In keeping with the choice for P(Sj^|Fc) above, 
one could choose P(Sj^|Tc) = P(Sj^) and P(Sj^|FQ) could be 
determined by the coverage characteristics of the detection 
system. In particular, P(Si|F^) = P(Fq | S • P(Si) /F^) where 
P(Fc) = 2 P(Fq| S j ) 'PCSj ) and the sum index j = 1,2, •• ,n. 

Using the correspondence between the events Fq and Dq n 
one can write P(FQ|Sj^) = 1 - (Pd)i* 

Both in the case of this procedure for negative information 
and the procedure for positive information, (Pd) i can be viewed 
as a measure of the specified performance of a detection system 
and both p and p can be viewed as measures of the degradation 
in its specified performance do to various operational factors. 

In particular, for p = 0, the expression for P(Sj^|C) is the 
same for the two negative information procedures if P(c|S£) is 
determined in the same way for both cases. In the view just 
expressed, this implies that the first procedure applies to a 
search in which a detection system is operated so that it 
achieves its specified performance. In the case of the positive 
information procedure, if p = 1, then the expression for P(S£|c) 
is the same as that which can be obtained from 

(33) P(Si|c) = P(c|Si) -P(Si)/P(C) 

where P(C) = 2 P(C | S j ) • P(S j ) and the sum index j = 1,2, ••• ,n 

if P(c|Si) is determined in the same way for both cases. 
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XV. Search Models and Search Theory 

Search theory provides a basis for determining optimal 
search plans for a target whose motion and location are 
determined within some bounds. Here, an optimal search plan is 
one for which the probability of finding a target within a given 
length of time is a maximum, the expected time to find a target 
is a minimum given the target is found or a search plan for which 
some other optimal search criterion is satisfied. 

Search theory results are based on models of the search 
process. To the degree that a search model describes a search 
process, an optimal search plan for a target that is based on the 
search model should provide guidance for the development of an 
operationally feasible search plan. However, because of the 
limitations of analytical search models, an optimal search plan 
that is based on an analytical search model may give only initial 
guidance in this regard. The optimal search plans that are 
described below illustrate this. The search plans are based on 
the random search model. Because of this, the requirement on the 
location of search track segments is not realizable and the time 
to resolve false alarms is ignored. 

Optimal search plans based on search models implemented 
through a monte carlo simulation are not considered here. 

However, with sufficient information, such plans have the 
potential of being both implementable and more optimal in a real 
sense than an optimal search plan based on an analytical search 
model . 
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Three Optimal Search Plans: The three optimal search plans 

differ through their definition of optimality. However, each one 
is based on the following search model; A target is fixed at 
some point in a region that consists of n subregions. A search 
in a subregion is a random search in the sense of the definition 
in Section X and a searchers sweep width there is a constant. 

In addition, a search of a subregion will not detect a target 
which is in another subregion. To determine a plan, let 
Sj^ = (the target is in subregion i} for i = 1,2, ••• ,n and 
let pj^ = P(Sj^) be the prior probability that the target is in 
the i^^ subregion. Let be the sweep width in the i^^ 

subregion. Let = Aj^/Wj^ where Aj^ is the area of the i^^ 

subregion and is the expected track length to find the 

target by a search of the i^^ subregion given the target is in 
the i^^ subregion, a characteristic length. The probability P 
that the target will be detected by a random search is given by; 



where the sum index i = 1,2, ••• ,n and Ij^ is the track 



The first optimality criterion is: Choose 1^ so that P 

is a maximum subject to the two constraints: 1. 1 = 2 Ij^ and 
2. 1^ > 0 where the index i=l,2, ••• ,n. Determining this 
choice is a nonlinear optimization problem whose solution is 
given in Reference 20. It is: 



(34) 



P = 2 [1 - exp(-li/6i) ] -Pi 



length of the search in the i^^ subregion. 



(35) 



lf/&i = ln(Pi/6i) - L(k) 

lt/8i = 0 



i = 1,2, 



i = k+1 , k+2 , • • • , n 



k 
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where L(k) = (1/S5 j ) *S [ 5 j • In (pj/6 j ) ] + 1/S5j and the sum index 
j = 1,2, ,k, where the subregions are relabeled so that the 

following order relation holds: Pi/Si > P2/^2 > *•* > Pn/^n 

where k is chosen so that for k+1 the solution for l^+i 
using L(k+1) is either negative or zero. 

The second optimality criterion is: Choose 1^ so that P 

is a maximum subject to the two constraints: 1. c = 2 Cj^ and 

2. cj^ > 0 where the index i = 1,2, ••• ,n, cj^ = k^-lj^ is the 

cost of the search in the i^^ subregion and k^ is the cost per 

unit track length in that subregion. For this criterion, the 
solution to the corresponding nonlinear optimization problem can 
be obtained from Equation 35 by replacing by £i = kj^*6j^ 

and labeling the subregions so that Pi/^i > P2/^i > *** > Pn/^i* 
The basis for this can be seen by replacing li/Si by its 
equivalent in the exponential term in Equation 32. 

The third optimality criterion is: Choose Ij^ so that the 

expected utility of finding the target is a maximum subject to 
the two constraints: 1. 1 = 2 Ij^ and 2. Ij^ > 0 where the 

index i = 1,2, ••• ,n. For this criterion, the solution to the 

corresponding nonlinear optimization problem can be obtained from 
Equation 31 by replacing p^ by q^ where qj^ = and Uj^ 

is the utility of finding the target given it is in the i^^ 
subregion. And, in addition, labeling the subregions so that 
<3l/Si > ^2/^2 > ••• > ^n/^n* basis for this can be seen by 

multiplying the summation term in Equation 34 by Uj^ so that 
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the resulting equation gives the expected utility of finding the 
target with the search. 

Equation 35 can be used to determine an order of search for 
the subregions which will effectively minimize the expected track 
length required to detect a target given it is detected. To do 
this, divide the available track length 1 into units small 
enough so that with a single unit only the subregion would be 

searched. Then allocate one unit to the search of the 1®'^ 
subregion. If the search is unsuccessful, determine the optimum 
allocation for two units. Then search with a second unit so that 
the first search with the first unit plus the second search with 
the second unit satisfy the optimum allocation for two units. If 
the search is unsuccessful, continue in this fashion until either 
the target is found or all the track length is expended. That 
this allocation order will effectively minimize the expected 
track length required to detect a target given it is detected can 
be argued as follows: Let L be the track length at detection, 

let 1^ be a unit of track length and let n be the number of 
units. Then the value of the probability P(L < i*lu) that the 
target will be detected on or before the i^^ step of the search 
for the given allocation order will be greater than or equal to 
its value for any other allocation order with the same allocation 
step size. Since the value of P(L < 1) will be equal to its 
value for any other allocation order of the optimum allocation 
and since P(L < i'l^lL < 1) = P(L < i*lu)/P(L < 1), the value of 
the distribution function FL(i*lu|L < 1) = P(L < i'l^lL < 1) 
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will be greater than or equal to its value for any other 
allocation order. This implies that the expected track length 
given detection E(l|l < 1) = S [1 - FL(i*lu|L < 1)] where the 
sum index i = 1,2, ••• ,n is effectively a minimum for the 
given allocation order. A search based on the optimum allocation 
given by Equation 35 and the given allocation order is 
equivalent to the following search: After an allocation of track 

length 1^ and an unsuccessful search, new values for P(Si) are 
calculated using Equation 29 and then Equation 35 is used with 
these new values to determine the next optimum allocation. A 
discussion of this procedure is given in Reference 6. And an 
example of its application is given in Reference 21. 

Equation 35 also defines an optimal search plan for a 
detection system that searches beams and can be described by 
Equation 33 by replacing Ij^ by t^ where t^ is the time the 
ith beam is searched and by replacing by where Tj^, a 

characteristic time, is the expected time to detect the target by 
a search of the i^^ beam given the target is in the i"^^ beam. 

For a more extensive discussion of search theory and its 
application to military operations research, see Reference 22. 
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